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Abstract. Bohmian mechanics allows us to understand quantum systems in the light of 
other quantum traits than the well-known ones (coherence, diffraction, interference, tunneling, 
discreteness, entanglement, etc.). Here the discussion focusses precisely on two of these 
interesting aspects, which arise when quantum mechanics is though within this theoretical 

■ framework: the non-crossing property, which allows for distinguishability without erasing 
^ ' interference patterns, and the possibility to define quantum probability tubes, along which the 

probability remains constant all the way. Furthermore, taking into account this hydrodynamic- 
like description as a link, it is also shown how this knowledge (concepts and ideas) can be 
straightforwardly transferred to other fields of physics (for example, the transmission of light 
along waveguides). 
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^ . 1. Introduction 

I At present, there is no doubt that quantum mechanics can be considered the most successful 

■ theory ever devised to explain the physical world. Its applicability ranges from very fundamental 
physical problems to the high-technology applications that are nowadays an important part 
of our daily life. This theory, though, still constitutes a veiled mystery at a deeper level 

^ ■ of understanding, for there is a lack of a clear interpretation of the physics underlying 

5^ i quantum systems. This is somehow connected to its widely accepted interpretation, namely the 

Copenhagen interpretation [1], which not only does not allow us to think of quantum systems 
as we do of classical ones, but it just forbids such a thing. 

A feasible way to surmount this drawback (although surely not the final one) comes through 
Bohmian mechanics [2-5]. According to this approach to the quantum world, systems are 
described in terms of trajectories evolving in configuration space. Rather than constituting a 
step backwards towards the classical paradigm, this quantum-mechanical formulation puts at 
the same descriptive level both classical and quantum systems, though the latter's behavior 
differs from the former's one through the dynamical particularities induced by the quantum 
laws of motion. For instance, these laws give rise to motions constrained to the configuration 
space, uncertainty principles, interference phenomena, discreteness, and other features which 
are absent in classical systems. Of course, in the statistical limit, i.e., when an appropriate 
sampling over many of such Bohmian trajectories is carried out, the results of the standard 
quantum mechanics are recovered. 

On the other hand, we find that in principle there are no restrictions to export the core 
idea in Bohmian mechanics, namely the possibility to reformulate quantum mechanics as a 
hydrodynamic-like theory, as well as the ideas emerging from it to other physical contexts where 
waves or distributions are the primary descriptor. This is somehow the other way around that 



Madelung [6] considered in 1926, when he tried to provide a hydrodynamical interpretation of 
quantum processes based on the analogy between the latter and classical hydrodynamical flows. 
Within this approach, Bohmian trajectories were just considered to be the "tracks" displayed 
by some hypothetical tracer particles left on top of the quantum fluid [7], something very similar 
to what can be found in modern studies of flows, where the trajectories displayed by a number 
of tracer particles is recorded and subsequently analyzed in order to extract information about 
the flow. By means of this method, for example, studies revealed anomalous diffusion and Levy 
flights in two-dimensional rotating flows [8, 9]. However, closer in spirit to the present discussion, 
we find the recent observations by Couder and coworkers [10-14] as well as by Bush [15], where 
a seemingly classical analog of Bohmian motion is found for droplets bouncing on a vertically 
vibrating bath of the same fluid. 

The purpose of the present work is to pose an analysis and discussion of some interesting 
properties of quantum systems and their consequences, which arise when they are investigated 
from a Bohmian perspective. Furthermore, the link that can be established between the Bohmian 
viewpoint and its application or extension to other non-quantum wave theories will also be 
discussed. In order to make clear the arguments exposed, the work has been organized as follows. 
In Section 2, the essential elements/concepts of Bohmian mechanics are introduced as well as 
some important consequences implied by this theory. The transfer of knowledge from Bohmian 
mechanics to other classical wave theories, such as electromagnetism, is discussed in Section 3 
in the context of light waveguiding. To conclude, a series of final remarks are highlighted in 
Section 4. 



2. What can be learnt from Bohmian mechanics? 

2.1. Some elementary background 

In the standard Bohmian approach [2], one usually starts considering the wave function in polar 
form, 

*(r,t) = i^(r,^)e^^("•'*)/^ (1) 

which is substituted into the (single-particle) non-relativistic, time-dependent Schrodinger 
equation, 

>n'4-(-l^'^vW (2) 



dt \ 2m 
to yield the system of real coupled equations 



rn 



+ V-[R' =0, (3) 



dt 

dt 2m 2m R ^ ^ 

In brief, the continuity (or conservation) equation (3) rules the ensemble dynamics, i.e., the 
number of particles described by the probability density has to remain constant; the quantum 
Hamilton- Jacobi equation (4), on the other hand, describes the time-evolution of the phase field, 
which governs the motion of quantum particles through the guidance equation 

vb = r = . (5) 

m 

In this way, the simple transformation described by (1) allows us to go a step beyond standard 
quantum mechanics. In the latter, ^ is all what is needed to describe and interpret the 
experimental outcomes; in Bohmian mechanics, on the other hand, the evolution of a tracer 



particle [7] according to the information supplied by ^ (the guiding field) is also considered to 
further understand the quantum system dynamics. 

Now, the same description can be formally inferred by simply assuming that the evolution 
of any quantum system can be identified with the diffusion of a (quantum) fiuid throughout 
configuration space, as formerly proposed by Madelung [6] . Therefore, as in the case of classical 
fluids, the elementary quantities of this approach would be the (probability) density of the 
quantum fluid and its associated probability current density, J, both related through the 
conservation (or continuity) equation 

| + V.J = 0, (6) 

where J = pwn and vj^ is a local hydrodynamic velocity field. The probability density and the 
probability current density are related to the wave function (1) by the transformation relations 

p = R\ J = i?2— , (7) 

m 

with the local velocity field being 

v„^i^^. (8, 

p m 

which is equivalent to the Bohmian one, v^, given by (5) (hence, from now on, both velocity 
fields will be denoted by v). The integration in time of (8) generates a family of streamlines 
or paths (for each wave function considered) along which the quantum fluid propagates, just as 
in the case of a classical fluid. As it can be inferred, the first equality goes beyond Bohmian 
mechanics and allows to define streamlines in any system characterized by a certain density and 
a vector that transports it through the corresponding configuration space, regardless whether 
such a density describes a quantum system or not. 



2.2. The non-crossing rule 

One of the most relevant properties revealed by Bohmian mechanics is that quantum fluxes 
cannot cross in configuration space. That is, at a certain time t, each point on configuration 
space has a uniquely-defined value of the local velocity field, except for regions where the wave 
function displays a node (i.e., ^' = 0). Therefore, two or more Bohmian trajectories cannot 
cross through such a point at the same time t or, equivalently, quantum systems cannot reach 
the same configuration at a given time following different quantum streamlines. In terms of the 
standard quantum mechanics, this just means that paths can be labeled, although they cannot be 
experimentally observed because this implies changing the process and therefore generating new, 
alternative paths. Obviously, this changes dramatically our physical view of the superposition 
principle, which states that any linear combination of solutions of Schrodinger's equation is also a 
formally valid solution to the same equation. Accordingly, given a coherent superposition of two 
counter-propagating wave packets, they can be evolved in time separately and then recombined 
at any time to determine when and where an interference pattern. Formally, there is no problem 
in proceeding in this way. However, from a (hydro) dynamical viewpoint the evolution of both 
wave packets has to be accounted for together for the interference patter to be observable 
(which is what actually happens in an experiment), this implying that interpretations of physical 
phenomena cannot directly rely upon the superposition principle. Bohmian mechanics provides 
us an answer for this based precisely on the non-crossing rule. 

In order to make clear such a statement, consider the wave packets above are represented by 
two free expanding Gaussian wave packets. 



(9) 



where the time-dependence of their spreading is given by 




crt = \^t\ = cro\ I + 2 • (10) 



The wave packet (9) propagates along the classical trajectory Xd = xq + vgi (with vq = po/m), 
while the associated Bohmian trajectories are given by 

x{t) = xa + — [x{0)-xo], (11) 

with x(0) denoting the initial position of the corresponding Bohmian trajectory. As can be 
noticed, there are two contributions to the evolution of the tracer particles, one coming from 
the classical drift and another one from a purely quantum origin [16, 17]. The former makes the 
particle to evolve according to only the laws of classical mechanics; the latter makes it to deviate 
in such a way that its motion accommodates to the evolution of the quantum fluid. The time- 
scale that separates both types of motion (to some extent) is given in this particular case by the 
characteristic time r = 2maQ/h [16, 17]. Depending on the relative value of t when compared 
with T, three regimes can then be found (in analogy with the regimes of optical physics [7]): 

• Ehrenfest or Huygens (r = 2m(jQ/h ;g> t): Bohmian trajectories follow classical ones (in 
agreement with Ehrcnfcst's theorem). 

• Near-field or Fresnel (r = 2maQ/h > t): Bohmian trajectories start to diverge from their 
classical counterparts. 

• Far-field or Franhofer (r = 2maQ/h <C t): Bohmian trajectories display a sort of rectilinear, 
uniform evolution due to the stationarity of the wave function far from any interaction 
potential. 

Therefore, by inspecting the topology displayed by the quantum trajectories, one readily notices 
that the quantum flow evolves from an initially confined fluid to a linearly expanding one, 
undergoing at times of the order of r a sort of internal boosting which bursts it open. However, 
this does not provide any clue on the true dynamics. 

In order to understand the importance of not disentangling the joint evolution of both wave 
packets, as indicated by the superposition principle, consider the full wave function, 

^'(r,t) = ^i(r,i) + V'2(r,t), (12) 

where the = p^^^e*"^'/^ are properly normalized and have the shape indicated by (9), with 
X2,o = —xifi = xo > and vi^q = — U2,o = fo > 0. It is straightforward to notice that 



p = Pi+ P2 + 2^pxpi cos (/p, (13) 
J = ^ [/9iV5i + P2SS2 + VpIpiV(Si + S2) cos + n (pl^Vps^^ - pT^pT^ sin v?] ,(14) 

with ip = {S2 — S\)/h. Substituting these expressions into (5) lead to the equation of motion 



/ l/2„ 1/2 l/2„ l/2\ . 

. _ 1 PiVSx + P2VS2 + VpIp^V(^i + S2) cos if h [Pi ^P-i -P2 VPi Jsmy 
m Pi+ p2 + 2y/p\p2 sin m pi + P2 + 2^ p\p2 cos (/5 

The crossed term in (13) and (14) arises from the interference of the probability densities and 
current densities associated with ^1 and ^2 after assuming the coherent superposition (12). 
However, as seen from (15), the local velocity field keeps a more complicated form, which cannot 




Figure 1. Bohmian trajectories associated with a 
coherent superposition of two Gaussian wave packets 
(gray hues) and with a single Gaussian wave packet 
(red hues). 
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Figure 2. Contour-plots representing the evolution in time of the probability density (a), the 
phase field (b) and the local velocity field (c) of a initial coherent superposition of two Gaussian 
wave packets. The Bohmian trajectories (black lines) associated with this problem are also 
displayed superimposed to each graph. 



be expressed as a simple superposition, thus indicating that the dynamics is going to be far more 
complex. This is illustrated in Fig. 1, where the single wave-packet Bohmian trajectories (red) 
diverge from the behavior displayed by the two wave-packet ones as they approach interference 
region, ending up eventually in the wrong place. This can be readily understood by looking at 
the three graphs displayed in Fig. 2: although the probability density (part (a)) does not provide 
any clue, both the phase field (part (b)) and the local velocity field (part (c)) immediately make 
apparent the dynamical consequences of the superposition. 

It is worth stressing that the non-crossing takes place regardless of the presence of the 
quantum potential. Rather it is an intrinsic or inherent property associated with the dynamics 
generated through a guidance condition. In classical mechanics we can also observe something 
similar, but in phase space, where phase-space orbits cannot cross through the same point at 
the same time. In such a case, the "guiding field" is the classical action (this was actually the 
idea that led Schrodinger in the derivation of his wave theory of quanta [18]). Nevertheless, in 
the two wave-packet problem presented above, the non-crossing can still be well represented by 
the action of an effective potential acting on a single wave packet, with the form of a repulsive 
wall and an attractive region, both depending on time [17]. 

2.3. Quantum probability tubes 

The possibility to study quantum processes in terms of Bohmian trajectories brings in another 
interesting property, namely the possibility to define quantum probability tubes [19], i.e., tubes 
in configuration space along which the integral of the probability density at a time t remains 
constant regardless of changes in their profile or in the shape of the enclosed sections of p. In 
order to show this, consider the so-called restricted probability, defined as 



Vit) = [ pir,t)dr, 
Jn 



(16) 



which gives the probabihty to find the system inside a certain region of interest of the 
corresponding configuration space at a time t [20-22]. From standard quantum mechanics (and 
making use of the divergence or Gauss-Ostrogradsky theorem), the variation with time of V{t) 
can be expressed as 

dV{t) f dp 



dt 



f %dv = - f (V-J)dr = - / J-dS, (17) 



where S is the boundary of $7. That is, the losses/gains of V{t) inside Q arc described by the 
outgoing/ingoing probabihty fiow through S, which is accounted for by the quantum current 
density J. Defining an arbitrary region ft in standard quantum mechanics is rather simple. 
However, how can this region be propagated in time so that V{t) keeps the same value along 
time? 

Although the answer to the question posed above can be quickly answered in classical 
mechanics by using arguments based on the Liouvillian structure of this theory [23], the 
same cannot be found in standard quantum mechanics, but in Bohmian mechanics and, more 
specifically, in its non-crossing property. As shown elsewhere [19], if the surface S is the geometric 
place formed by a (discrete or continuous) set of initial conditions, because of the non-crossing 
property, any probability enclosed inside such a set will remain constant at any subsequent 
time. This is based on the fact that at any given time t no trajectory inside or outside the 
boundary defined by the set of separatrix trajectories will be able to penetrate into or exit from 

respectively. Therefore, if the probability P{t) is measured in terms of the accumulation of 
trajectories embraced by E, i.e., 

N 
i=l 

this number must remain constant at any time provided S(i) is a causal map of S(0) at time 
t. An important consequence of this result is that any restricted probability can be determined 
directly from the initial state if the end points of the associated separatrix trajectories are known. 
This means that, in principle, one could determine (or, at least, estimate) final probabilities 
without carrying out the full calculation [24], but directly from the particular region covered 
by the initial wave function causally connected with the feature of interest (e.g., a transmitted 
amplitude or a diffraction peak), since 

Voo= [ p{roo)dVoo = hm / p[Y{t)]dY{t) = lim V{t) = [ p[ro]dro = Vq. (19) 

Here, p(roo)iiroo is the probability to find the system confined within the layer of configuration 
space defined by the separatrix trajectories ending at r^o and those ending at = [r -|- (ir]oo- 
Analogously, p[Y{t)]dr{t) is the probability to find the system confined within the layer defined 
by the sets of separatrix trajectories r(i) and r'(t) = [r + cir](i) at a time t. Note that, since 
the probability density is being evaluated along quantum trajectories, it is preferably denoted 
as /9[r(t)] instead of p{r,t), with Too = r{t — > oo). 

In order to illustrate this result, consider the scattering of a wave packet with a rather 
arbitrary shape by a barrier, as displayed in Fig. 3(a) (details on this calculation can be found 
in [19]). The initial and final probability densities are denoted by the black and red solid lines, 
respectively, while the barrier is indicated by the blue shadowed region. The propagation goes 
on until the probability in the region covered by the barrier (intra-barrier resonance) becomes 
negligible, which can be better seen in Fig. 3(b), where the transmission (green dashed line), 
reflection (blue dash-dotted line) and intra-barrier resonance (red dotted line) probabilities are 
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Figure 3. Scattering of a general wave packet by an almost square potential barrier, (a) 
Initial (black line) and final (red line) probability densities; the barrier is denoted with the 
blue shadowed region, (b) Time-dependence of different restricted probabilities: transmission 
(green dashed line), reflection (blue dash-dotted line) and intra-barrier resonance (red dotted 
line) probabilities. The black solid line denotes the transmission probability directly computed 
from the initial state (see text for details), (c) Bohmian trajectories illustrating the tunneling 
dynamics. The separatrix trajectory is displayed with the red thicker line, (d) Initial probability 
density split according to the separatrix initial condition: initial conditions in the green shadowed 
region contributes to transmission and those from the blue shadowed one to reflection. 



recorded along time. These are the three restricted probabilities considered in this example, 
with Jltrans being the region beyond the right edge of the barrier, Qj-es the region bound between 
the two edges of the barrier, and r^rcfl the region to the left edge of the barrier. After t ~ 1.15, 
"Pj-gs PS and Ptrans reaches its maximum, asymptotic value; Vrefi, on the contrary, reaches its 
minimum earlier, at t ~ 0.75. This means than T^oo, trans can be computed by simply choosing 
as the separatrix trajectory a Bohmian trajectory that reaches the left tail of the transmitted 
wave packet in part (a) and propagating it backwards in time, as seen in part (c). Integrating 
the probability density between the position of this trajectory at each time and infinity (for 
practical purposes), it is then found that Poo.trans remains constant all the way, as the black 
straight line shows in part (b). Regarding the initial set (probability density), it is found (see 
part (d)) that any initial condition started to the right of the separatrix initial condition will 
get into the transmitted region, while if it starts to its left, the trajectory becomes reflected. 

3. Where can this knowledge be transferred to? 

It is very common to separate quantum physics from other wave theories based on the difference 
in the treatment of probability amplitudes with respect to other types of waves (which are 
considered as perturbations propagating through a given medium or vacuum). However, the 
way how these wave theories operate is in essence the same, for both include the same type of 
concepts (e.g., coherence, interference, diffraction, tunneling, discreteness, etc.) and principles 
(e.g., uncertainty, superposition). Taking this into account, it is possible to establish a sort 
of feedback regarding the understanding of both wave phenomena, in general, and quantum 
phenomena, in particular, through the hydrodynamic (Bohmian) approach considered above. 



To briefly show how this transfer of knowledge can be done, consider the case of 
electromagnetic fields, for example. Bohmian mechanics was formerly formulated to describe 
massive particles, resulting very appealing to explain event-to-cvcnt experiments like those 
carried by Merli et al [25], Tonomura et al. [26], or Shimizu et al [27]. Now, there are also 
event-to-event experiments with light, such as those performed by Dimitrova and Weis [28], 
which claim for a treatment on equal footing. As shown by Prosser [29] , this is actually possible 
by directly considering Maxwell's equations (which put electromagnetism at the same level 
of Schrodinger's wave mechanics [30]). In this case, given an electromagnetic field^ defined 
by an electric field E(r) and a magnetic field H(r), the two key elements necessary to define 
the corresponding streamlines are the electromagnetic energy density and the Poynting vector 
(electromagnetic current density), i.e., 

C/(r) = \ [6oE(r) • E*(r) + MoH(r) • H*(r)] . (20) 
S(r) = iRe[E(r)xH*(r)], (21) 

respectively. Since the electromagnetic energy density is transported through space in the form 
of the Poynting vector, a local velocity field can be defined [31] in analogy to (8), which reads 
as 

S(r) = C/(r)v, (22) 

from which electromagnetic energy fiow lines or photon paths are obtained by integrating the 
equation 

* - 

ds~ cu{vy ^ 

along the arc- length coordinate s (which can be referred to a proper time r = s/c, with c being 
the speed of light). By means of this procedure, the interference of diffracted polarized beams 
was studied [32,33], showing a good agreement with experiments carried out later on to infer 
the average paths displayed by photons [34]. 

The previous discussion focussed on Maxwell's equations. However, it is not necessary to 
stay at such a level in order to define hydrodynamic streamlines. As mentioned above, this 
can be done at any level where the elementary descriptive tool is a wave equation. Thus, 
consider, for example, the transmission of a light pulse through a waveguide within the small- 
angle or paraxial approximation [31]. Assuming that the optical axis is oriented along the z-axis 
and the electromagnetic field passing through the waveguide is time-harmonic, the field can be 
approximated by a plane wave along the z-direction modulated by a certain complex-valued 
amplitude, 

$(r) = (^(r)e^'=^^ (24) 

where kz = nok, k = 27r/A, A is the light wavelength in vacuum and no is the bulk refractive 
index. Substituting this expression into Helmholtz's equation, 

+ n^k'^^ = 0, (25) 

with n being the refractive index inside the waveguide, one readily finds 

2ikz ^ + ^ = - Vi<^ + ik! - n'k')<f>, (26) 



^ Timc-indopcridcrit electromagnetic fields will be considered here for simplicity, although time-dependent ones 
should be used, in general. 
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Figure 4. Light transport through a waveguide with the shape of a Y-junction. (a) Initial 
(black) and final (red) probability densities, (b) Optical streamlines illustrating the dynamics 
inside the waveguide, (c) Enlargement of part (b). 



where is the transverse part of the Laplacian (for example, in Cartesian coordinates, it reads 
as V^^ = d'^ jdx^ + d"^ /dy^). If, apart from the paraxial approximation, it is also assumed that 
the slowly varying envelope approximation holds (i.e., the amplitude (p varies slowly in space 
compared to 2-K/kz), the highest-order derivative in z can be neglected (i.e., d'^(\>ldz^ w 0) and 
Eq. (26) can be recast as 

This equation is, indeed, isomorphic to Schrodinger's one, with z playing the role of the 
evolution parameter (rather than the time t) and = n^h the role of an "optical mass". 
Equation (27) has been used, for example, to study the design of waveguides with optimal 
conditions of light transmission [35-37]. However, regarding our discussion here, it also brings 
the possibility to define a family of optical streamlines if the pulse (j){r) is expressed in polar from, 
as ^(R, z) = p^/^(R, z)e*'^(^'^\ where R denotes the vector perpendicular to the propagation 
direction (z). In this case, the local velocity field reads as 

which gives rise to such streamlines when integration over z is carried out. 

In Fig. 4 we have an example of how these trajectories look like in the particular case 
of a Y-junction [38], which splits the initial pulse (red line) into two emerging pulses (black 
line), as seen in part (a). As it can be noticed in part (b), also here the electromagnetic flux 
follows the non-crossing rule, which allows to distinguish between different paths (and therefore 
to also define density tubes) without disturbing the experiment, also in agreement with the 
experiment performed by Kocsis et al [34]. Actually, in the enlargement of part (c), a behavior 
mimicking that of the frustrated internal reflection (inside the waveguide) can be observed, 
although contrary to geometric rays, the optical streamlines follow the wavy behavior dictated 
by the pulse inside the wave, thus being affected by diffractions and interferences. 



4. Final remarks 

We are used to consider that quantum mechanics is about small objects. However, this is just a 
misunderstanding, for quantum mechanical behaviors arise whenever objects are simple enough, 
i.e., whenever the action of any internal or external surrounding degrees of freedom does not 
play any relevant role in their dynamics. As an example, just consider the interference patterns 
that can be obtained from large organic molecules [39] or Bose-Einstein condensates [40,41]. 
Such a misunderstanding comes from technological lacks in the past, when quantum mechanics 



was developing, as well as black-box interpretations, like the Copenhagen one. The possibility to 
carry out interference diffraction of large objects, particle by particle, requires of new viewpoints 
and interpretations that may go beyond the blind ones available at present. 

In this sense, the role of Bohmian mechanics as a theory of quantum motion that describes 
the evolution of tracer particles when traveling along quantum fluids (i.e., the fluid associated 
with the diffusion of the wave function in a configuration space) are highly needed. Even if 
at a statistical level predictions are the same as those already provided by standard quantum 
mechanics, the fact that hydrodynamic-like analyses can be carried out, elucidating properties 
such as the flux non-crossing, results very appealing, for it allows us to understand the (quantum) 
physical world in a different manner, out of conceptual restrictions. 

Furthermore, such approaches should also be extended to other fields of physics, for they could 
help us to get different insights of already well-known and established processes and phenomena. 
This is the case, as it has been seen, of the so-called classical electromagnetism. In this case, 
not only theoretical descriptions have been developed to better understand the evolution of 
the electromagnetic energy through space, but also experiments have been motivated that have 
corroborated the models proposed [34]. Of course, and there is no doubt about the important 
role these approaches present from a pedagogical viewpoint [42] by allowing us to connect the 
distribution of energy (probability) through space with its diffusion (current densities). 

Finally, I would like to close this work with a simple, but striking reflection. Usually, one is 
always tempted to think that only those physical problems that are directly connected with a 
particular practical application are of interest and, therefore, worth spending time and efforts 
to investigate them. Of course, this is a legitimate point of view regarding where efforts and 
resources should be targeted to. However, this strategy cannot be used to get deeper into 
the physical theories and approaches behind such problems in order to getter a better and 
deeper understanding of them. This is precisely what happens with quantum mechanics and, by 
extension, with wave theories. Nowadays nobody questions their use, but there are less and less 
researchers facing the problem of the understanding of the quantum world and the paradoxes it 
brings in to our intuition, which in spite of all is shaped by a classical world. 
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